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Using the framework of the Chromo Dielectric Model we perform an analysis of color electric 
flux tubes in meson-like qq and baryon-like qqq quark configurations. We discuss the Abelian 
color structure of the model and point out a symmetry in color space as a remnant of the SU(3) 
symmetry of QCD. The generic features of the model are discussed by varying the model parameters. 
We fix these parameters by reproducing the string tension r = 980MeV/fm and the transverse 
. . . width p = 0.35 fm of the qq flux tube obtained in lattice calculations. We use a bag constant 

' B^/* = (240 - 260) MeV, a glueball mass nig = (1000 - 1700) MeV and a strong coupling constant 

\ Cpas = 0.2 — 0.3. We show that the asymptotic string profile of an infinitely long fiux tube is 

■ already reached for qq separations R > 1.0 fm. A connection to the Dual Color Superconductor 
(N ■ is made by extracting a magnetic current from the model equations and a qualitative agreement 

^ I between the two descriptions of confinement is shown. In the study of the qqq system we observe 

Q , a A-like geometry for the color electric fields and a Y-like geometry in the scalar fields both in the 

' energy density distribution and in the corresponding potentials. The resulting total qqq potential is 

' described neither by the A-picture nor by the Y-picture alone. 
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■ I. INTRODUCTION 

o : 

, The structure of hadrons is still a subject of discussion. It is widely accepted that Quantum Chromodynaniics 
(QCD) is the right theory for strong interactions, and that the properties of hadrons can be described within this 
i-C framework. QCD has been tested successfully in the region of large momentum transfer, where perturbative methods 
work due to asymptotic freedom. However, in the region of small relative momentum, where the formation of hadrons 
O i' sets in and confinement plays a dominant role, a calculation from first principles is still limited to lattice techniques. 
O . In order to understand the formation of hadrons out of quarks and gluons dynamically, and in turn the interactions 
' of hadrons with each other, one still has to rely on models that include the phenomenon of confinement. Such 
^ , phenomenological models are for example the bag model [mB, where quarks are treated as free particles restricted 
k>( ' predefined bags, the quark molecular dynamics model y,y|, where colored quarks obey the classical Hamilton 

rS ] dynamics bound by a linear confining potential, or the model of the dual superconductor |^|^, where confinement 
is achieved by monopole condensation and an accompanying dual supercurrent. The model of the stochastic 

vacuum relies on the calculation of Wilson- loops in a Gaussian approximation 9, 10] which leads to a linear potential 
between quarks and anti-quarks. 

In this work we adopt the Chromo Dielectric Model (CDM) [Hill 13, which is an extension of the bag model in 
the sense that bags are formed dynamically in the presence of quarks. Mack j3| and Pirner et al. 0, Ha . IT^ have 
given a renormalization group derivation of the lattice colordielectric model, which also has a scalar field modelling 
confinement, but keeps strongly coupled non-Abelian fields in the large distance action. A Monte-Carlo calculation 
within the lattice colordielectric model was done in [isl Hgj . The model has already been used to calculate hadron 
properties like low lying baryon masses, the nucleon magnetic moments and the (axial- vector) /vector coupling constant 
ratio 20] and nucleon-nucleon interactions in vacuum and in nuclear matter 15,21]. 22, 23, 24]. In another approach 
the description has been used within a transport theoretical scheme to describe the dynamics psl [26l l27l | of quarks 
bound in nucleons and strings. In l28l] a full molecular dynamics simulation for colored quarks was performed, showing 
the ability of the model to produce color neutral hadrons out of a gas of colored quarks, thus giving for the first time 
a microscopic description of hadronization from a quark gluon plasma. 
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The parameters used in |25|, l26l 1271 l28l |. which define the model, were mostly motivated by phenomenological 
arguments and not subject to a further investigation. The resulting color flux tubes are rather large objects with a 
string radius up to 1 fm. In addition the linear rising qq potential was only seen for quark separations R > 1 fm. On 
the other hand, lattice calculations 29, "sS] indicate, that the radius of a colored flux tube is much smaller than 1 fm, 
and that the qq potential already develops a linear rising term for separations larger than R = 0.2 fm. In addition, on 
the lattice a clear Coulomb-like potential was observed for quark separations R < 0.2 fm, which has not been resolved 
in '28l. In the present work we model the results from lattice calculations as well as possible within the framework 
of the CDM. On the lattice the most accurate results were obtained in SU(2) We assume that the shape of the 
color flux tube does not depend significantly on the underlying Lie algebra. Therefore we compare the results of our 
calculations to those obtained in lattice SU(2) theory. The quantities that we want to reproduce are the transverse 
sha pe o f a flux tube of given length, and the linear coefficient of the qq potential obtained in meson spectroscopy 

Having established a set of parameters matching the criteria above, we then analyze the structure of the flux tubes. 
The emphasis will be on the question how the string is build up when the constituents of the string are separated 
from each other. In varying the qq distance R we probe both the perturbative (small R) and the non-perturbative 
region (large R). We will also see how fast the transition from one to the other sets in and from which distances on 
the string picture holds. 

The formation of strings is not restricted to qq objects. In QCD one has the possibility to build up color neutral 
objects from three quarks. When the pairwise quark separations are large compared to the characteristic width of the 
string, color flux tubes will stretch between the quarks. In general two geometrically different pictures are possible. 
The first is the so called Y-geometry, where three fiux tubes meet at a central point [33. Issj . The second one is 
the A-geometry, where the three quarks are connected pairwise ISo. [371 . The three-quark potential emerging from 
these two pictures has been compared to lattice results in [s^ and |39l |. Due to the lack of numerical precision the 
two groups obtained different results. Therefore lattice calculations cannot clearly rule out one of the pictures so 
far. Within our model we are able to describe those baryonic quark configurations not only on the level of the qqq 
potential but also on the level of the energy distributions. We can study the structure of the formed flux tubes and 
discriminate between the two geometries. 

The main goal of this work is to fix the model parameters on lattice data of qq flux tubes. The parameters obtained 
will be used to describe both the shape and the potential of the three quark system within CDM. The structure of 
this work is as follows. In section^ we present the model and its specific mechanism of confinement. In section Hill 
we analyze the dependence of the shape and the potential of a qq string on the variation of the parameters introduced 
in the model. We compare our numerical qq results for three specific sets of parameters to those obtained within 
lattice gauge calculations in section Hvl and extend the analysis to baryon-like three-quark systems in section In 
the appendix we give a description of the algorithm used in our numerical calculations. 

II. THE CHROMO DIELECTRIC MODEL 

Presumably the non-Abelian gluon interactions of QCD are responsible for a highly structured non-perturbative 
vacuum. Although it is difficult to disentangle these interactions from first principles, the large scale behavior of 
strong interactions might be simple. In the Chromo Dielectric Model (CDM) one assumes, that the QCD vacuum 
behaves as a perfect dielectric medium, i.e. as a medium with vanishing dielectric constant. Colored quarks embedded 
in this vacuum produce electric color fields. In the presence of the dielectric medium these fiels are compressed in 
well defined flux tubes connecting quarks with opposite charge. 

To be more specific, the medium is described by a colorless scalar field a which mediates the vacuum properties via 
the dielectric function k((t). The confinement field itself is evolving in the presence of a scalar self interaction U{a). 

As the dynamics and the non-Abelian interactions of the gluon sector are merged in the confinement field a and 
its dielectric coupling K{a) one is left with a set of two Abelian gluon fields A'''" that interacts with the dielectric 
medium. In principle it is possible to formulate the model with dynamical quarks, described by a Dirac-like Lagrangian 
pol I23. liof . but in this work we concentrate on the structure of color flux tubes. Therefore we treat the quarks as 
external sources of color fields. The CDM can now be defined by the following Lagrange density: 

C = Cg + C„ , (la) 
Cg ^ ~lK{a)F^,F^^'''--gs3;,A'^'^ , (lb) 
= ^d^ad^a~U{a) , (Ic) 

pf.i^,a ^ Qt.j^i'M _Q,.j^f.,a^ a (={3, 8} . (Id) 

The color fields couple to the color charge current = (p"j j°) classically with a strong coupling constant gs 
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which we have not included in the definition of the current. The charge density is defined as a sum over all charged 
sources, i.e. the quarks, p°'{x) = q1w{x — Xk), where the are the color charges of the quarks. In the static case 
we are interested in, the spatial part of the color current vanishes j = 0. The quarks are in principle point-like objects 
but in our numerical analysis described in sec. lllll we assign a finite Gaussian width w{x) = (27rrQ)~^/^ exp(— a?^/2rQ) 
to each quark. The width is introduced for numerical reasons (see appendix lA 1|) . Throughout this work we adopt 
a value of ro — 0.02 fm, which is large enough to resolve the Gaussian distribution 'w{x) and small compared to the 
dimensions of the flux tubes. 

As we are working in an Abelian model inspired by QCD, we have three different colors interacting with only two 
Abelian color fields. In the color space we choose an arbitrary but fixed base for the quarks in the fundamental 
representation, |red) = (1,0,0)-^, |green) — (0,1,0)-^ and |blue) — (0,0,1)-^. The color charges q"" are then defined 
as the diagonal entries of the corresponding generators i° of the color group in the same representation g° — {c\t°-\c), 
where c S {red (r), green (g), blue (b)} and a S {3,8}. The numerical values of the color charges can be read off 
from table|l]and are depicted in figure^ The generators are normalized according to Tr f^t^ = 8°"^ jl. Note that by 
this definition the qa are reduced by a factor of two as compared to and that the strong coupling parameter gs is 
enhanced by the same factor leaving gsqa fixed. 



color 






red 


1/2 


l/(2^/3) 


green 


-1/2 


1/(2^3) 


blue 





-l/%/3 



TABLE I; The color charges g° of the three colors with respect to the two Abelian color fields. 



2^/3g' 




2g3 



FIG. 1: The color charge g^'* with respect to the color fields A?'^ . 

In the Abelian projected theory one is left with two independent color fields and a remaining U(l)xU(l) gauge 
symmetry. The corresponding gauge transformations are 

|c) ^ i7(2;)|c) =diag(e^9=xi^e^9.X2^g-«s,(xi+X2)j (2) 
^ {Al + d^Ba)ta, (3) 

where Xi = % + and X2 — —^■ + '2^' '^^^^ means that the color charges (f"!^ are conserved independently and 
in turn the color fields F^^''^ would be observable fields. However in the Abelian approximation there is a further 
symmetry, namely a symmetry under special discrete rotations in color space 

\c)~.\d)=V\c) . (4) 

In matrix form the color rotations are given explicitly as: 
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: 




= ei, 




= £2, 


^33 


= £l£2 


V2 : 




= ei, 


Vis 


= £2, 




= £l£2 


^3 : 


''is 


= ei, 


^21 


= £2, 


^^32 


= £l£2 


Vi : 


^11 


= ei, 


*^23 


= £2, 


^32 


= -£l£2 


^5 : 


^13 


= ei, 




= £2, 




= -£l£2 


^6 : 




ei, 


^21 


= £2, 


^33 


= -£l£2 



(5a) 
(5b) 
(5c) 
(5d) 
(5e) 
(5f) 

where €1/2 — ±1 and all other matrix elements being zero. There are therefore 4 different copies of each of the Vi 
differing from each other in the signs ei/2- These rotations V act either as a cyclic exchange of all colors r g b 
(j5b|) and b ^ g ^ r (|5c|l or as a pairwise exchange of two colors g ^ b Ij5d|l . r ^ b lfSe|) and r ^ g with additional 
phases £1/2- The transformations themselves form a global subgroup D of SU(3), but are not independent from the 
former U(l)xU(l) gauge group in the sense that out of the set of four copies of Vi three of them can be constructed 
by a combination ViU{x) with U{x) S U{1) x U{1). The discrete color rotations Vi transform the gauge fields A^-^ 
and A'^'^ only into each other without mixing to the non-Abelian gauge fields: 

= ^M,3^^ ^M,8^g ^ VA^l/t = ^'M,3^^ ^ A'^-'hs, (6) 

with (A^'3)2 + (^M.8)2 ^ (yi'M,3^2 ^ (^'^",8)2 ^/^.a ^ ^M^a. The samc is true for the color fields F^''-^ and F^''-^. 
They are therefore not invariant under the rotations V G D. This is a relict of the full SU(3) gauge symmetry. Of 
course, the action density |^ and the corresponding energy density, which are the only physical meaningful quantities 
in the model, are invariant under V G D. 

In the confinement part of the Lagrangian (|lc|l the scalar self interaction U{<7) is of a quartic form, i.e. 

U{a) = B + a(T^ + b(7^ + ca*. (7) 

The form is chosen to develop two stable points. A metastable one at cr = and a stable one at the vacuum 
expectation value a = tivac (see fig. |2J • The requirement that U has an absolute minimum at cr = cr vac , together with 
[^(fvac) = 0, leaves only two additional free parameters which we choose to be the bag constant B and the curvature 
TOg = ?7"(crvac) of the potential J7(cr) at the absolute minimum. Since the confinement field cr absorbs non-Abelian 
gluon properties of QCD, we can interpret rUg as the mass of the lowest collective gluon excitation, i.e. the glueball 
mass. The parameters a, b and c can therefore be expressed by the quantities B, nig and cr vac: 

1 mgaj, ~ 12B 

a = 7^— 2 , (8a) 

"l^O'vac - SB , ^ 

b = , 8b 

vac 

1 "^gf^vac - 6B 
^ "^vac 

In order to have a local minimum at cr = we must fulfill a > or 

mlol^, > 12B (9) 

The generic form of U (cr) is shown in fig. |21 For the equality in eq. the potential U has only an infiection point 
at cr = (solid curve). For fixed nig and cr vac we can vary B up to an upper limit given by eq. ^ (dashed curves). 
Alternatively we can fix B and cr vac and vary Wg starting from a lower bound given by eq. lO (dash-dotted curves). 

The two (meta-)stable points separate two different phases of the vacuum. In the perturbative vacuum at ct = 
where the dynamics is driven by short range interactions electric fields can propagate freely, i.e. the vacuum is 
described by a dielectric constant k — 1. In the non-perturbative phase cr ~ cr vac the dynamics is dominated by 
long range interactions and the dielectric constant nearly vanishes {k = Kvac ^1)- In the limit Kvac the 
non-perturbative vacuum behaves as a perfect dielectric. Between the two vacua the dielectric function k{(j) drops 
continuously from 1 to Kvac (fig- El- We choose for its parameterization a 5th order polynomial 

1 + hs^ + kis^ + hs^ , < s < 1 
<s)^{ 1 , s<0 , (10) 

'^vac 7 S > 1 
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FIG. 2: The scalar potential U{a) with mg fulfilling the equality in eq. @ (solid curve), fixed nig (dashed curves) and fixed B 
(dash-dotted curves). 

with s = cr/cr vac and with coefficients 

fcg = i(29Kvac-20) fcg = "10 

h = (15-23Kvac) ""^^ ^ 15 . (11) 

fc5 = i (19Kvac - 12) ^5 = - 6 

The coefficients are chosen such that the first two derivatives of k at cr = vanish and both k and its first two derivatives 
at (T = CTvac are proportional to Kvac and thus vanish in the Hmit /Cvac — > 0. The form of this parameterization is only 
weakly sensitive to the value of Kvac- We will see that physical quantities do not depend on the actual value of Kvac 
once it is chosen small enough. In fig. |31we have chosen a value k^^c = lO"** and the difference to Kvac = 10^^ is 
hidden within the linewidth. Of course other parameterizations are possible. In |4lj | we used K{a) = k~^^ but this 
functional form depends more strongly on Kvac and thus has more influence on physical quantities. The dielectric 
function approaches Kvac faster for decreasing Kvac- As a consequence the transverse shape of the color flux tube 
becomes steeper with vanishing Kvac and the energy of the system still increases even for very small values of Kvac- 
The actual polynomial choice is very similar to that in [23| and 0|- All parameters defined in eqs. Q), (0, ((HJ are 
subject to an investigation described in sec. lIIIBl The choice of Kvac is discussed in appendix lA 21 

From a variational principle we can derive the equations of motion for the gluon fields = (</>", A") and the 
scalar confinement field a: 

d^inF^^-'n = , (12a) 

d^d^-a = -U'{a) - y{a)F^,F^''^^ , (12b) 



where the prime denotes a differentiation with respect to a. The color field tensor F'^^'^ "' in eq. (|ld|) determines the 

color electric and magnetic fields Ef = -i^"^" ^ (-V^" - dtA'^'j and Bf = -^e^^juF^^'"- = (V x A'')^ respectively. 

With the help of the electric and magnetic fields we can recast eq. H12a|l into the two sets of inhomogeneous Maxwell 
equations, which we supplement with the two homogenous ones, that are fulfilled automatically by the definition of 
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FIG. 3: The dielectric function drops from its perturbative value k = 1 at cr = to its non-perturbative value Kvac ^ 1 at 
o" ~ o"vac. Here we have used Kvac ~ 10^*. 



the field tensor in ljld|l 



V X ff" - dtD" = gsj" 

V X i?" + 9(5" = 

V-B = , 



(13a) 
(13b) 
(13c) 
(13d) 



where we have introduced the electric displacement D'^ — K{a)E"' and the magnetic field iJ" — K{a)B°-. The energy 
of the system is given by 



(14) 



In this work we are interested in static solutions of given quark configurations, i.e. /= 0, and one can assume that 
j^fj,,a g^jjj ^ g^j-g g^j^gp ^jjjjg Indep eudcut . In this case the energy is minimized with B'^ = and eq. ()14|l reduces to 



Etot — 


Ed + -Evol + Esur 


(15a) 


Ed — 


i J (frE" ■ D'' 


(15b) 


Evol = 


j d^rU{a) 


(15c) 


Esur — 




(15d) 


subject to the constraint in form of Gauss's law (|13a|l 


The equations of motion now read 




V • (k(ct)V0") 




(16a) 


VV-C/'(a) 


= -^\ti'{a)E'' ■ E'' , 


(16b) 



Note that Gauss's law has to be fulfilled always and that it is conserved by the dynamics of the system. Note also 
that in the electric part of the energy we have included the self energy of a charge distribution which diverges for 
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point-like particles. The absence of the magnetic field might be due to the neglect of quantum effects. On the lattice 
[2^23 it is seen, that energy density and action density differ to some extent, i.e. (i?")^ ^ 0. As pointed out in 
in the framework of the stochastic vacuum model, the squared magnetic field is dependent on the renormalization 
scale. In the cited work the calculation of a qq flux tube was performed at a scale, where the magnetic fleld vanishes 
as well. However, in the work of Pirner in a non-Abelian version of the dielectric model |17| a non- vanishing magnetic 
field arises naturally. 

The coupled system of equations (|16|l for the two color electric potentials 0" and the confinement field a are the 
fundamental equations to be solved in this work. This is achieved numerically by the Full Approximation Storage 
algorithm described in appendix lA II 

There is one comment about the validity of the classical treatment of the model. By describing the quark fields 
by classical charge distributions one might expect some modification of the results due to the neglect of quantum 
mechanical aspects. One famous modification of string-like objects is for example the Liischer-term showing up in the 
gg-potential |4J| and corrections to it 45], which is due to quantum fiuctuations around the classical solution. We 
will discuss the gg-potential within our model later in scc. lIII Al Another modification might be expected due to the 
neglect of possible superpositions of states in color space. We will show now, that in the Abelian approximation the 
energy density is unaffected by changing from a quantum mechanical superposition to the classical analog used in our 
calculations. 

To be explicit regard a meson type state \qq) and a baryon type state \qqq) 



(17a) 



\qqq) = |B) = ^ e,jfe|(c,V'i) (cjV'2) (cfeV's)) , (17b) 

ijk 

where c(c) denotes the color of the (anti) quarks and '0(V') tbe spatial part of the (anti) quark wavefunction. Here 
we have neglected all other quark quantum numbers like spin and flavor as they do not enter the Lagrangian. The 
contraction with the total anti-symmetric tensor eijk in eq. Ijl7b|l ensures the anti-symmetry in the color part of the 
baryonic wavefunction. 

The classical analogs of these states in our model are 

\qq) - |A/ei) = |(r^)(f^)) (18a) 
\qqq) = = |(r^i) (5^2) (6^-3)) (18b) 

where r, (f), g, b are the (anti-)colors of the particles. The explicit choice of colors indeed is irrelevant due to the global 
color symmetry discussed before. The particles are supposed to be located at different positions, i.e. |M(A/ci)) and 
\B{Bci)) describe extended objects. 

The charge density p°'{x) is the expectation value of p^{x) = p„{x)t'^ in a given hadronic state, where Pn{x) — 
Si^x — Xn) is the one-particle density operator. The sum runs over all quarks and anti-quarks and the index n indicates 
that the operator acts on the nth particle in the state. 

It turns out that the charge densities of the quantum-mechanical states H17|l vanish everywhere, whereas they are 
finite for the classical analogs. Note that the total color charge vanishes in both cases. However, the charge density 
is not an observable quantity due to the global color symmetry and one better should look at the energy density. If 
one rewrites the electric part of the energy (|15bp in terms of the charge density one gets: 

E., = U d^xdV^^P^. (19) 
J \x~v\ 

Here we have used the perturbative expression (k — 1) and (/)°(a;) — j d^yp°'{y) /\x — y\. The modification of the 
electric energy due to the scalar function «;(a) in the non-perturbative case does not change our argument. Thus the 
interesting part of the electric energy is given by the expectation value of the squared charge operator p"" {x)'p°' (y) . In 
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the Abelian approximation (a € {3, 8}) it turns out that 



+ (^|p(2^)p(y)l^>] 

= (Mei|p«(x)p«(y)|Mei) (20a) 

3 

{BirixWivm = ^(^„ip(x)p(y)i^„) 



ra=l 
3 



2 

{m,„} = l 

= {B,i\r{x)r{y)\Bci) ■ (20b) 

Here C^'' = 1/3 is the eigenvalue of the quadratic Casimir operator in the fundamental (3-dimensional) representation 
in the Abelian approximation defined by X)a ~ C'f^Is? ^ G {3,8}. The first term in both equations denotes the 
self energy of the particles and the second one the two particle interaction. In the baryon case, the interaction is 
accompanied by an additional color factor 1/2 due to the interaction between two quarks whereas in the meson case 
it is a quark-antiquark interaction. Note that the equalities in eqs. ()20ll between the quantum-mechanical and the 
classical expressions are only valid in the Abelian approximation where a G {3, 8}. The interaction in the states \M) 
and \B) is accompanied by a color factor — Tr f^f^/N, where N is the normalization constant appearing in eq. H17|l. 
In the classical analogs the same factor amounts to —{c\t°'\c){c'\t°'\c'). By explicit calculations one sees, that these 
expressions are the same when summing over a e {3, 8} but differ by a factor of four when summing over a £ {1 ... 8}. 
We conclude that the classical treatment of qq and qqq states is reasonable in the Abelian approximation as it does 
not influence the observable energy density. 

III. ANALYSIS OF THE MODEL 

The confining properties of the model are ruled by the interaction of the color fields with the dielectric medium 
psLli^ . In the absence of any colored quarks, all color fields vanish and the confinement field will take on its vacuum 
expectation value cr vac- If quarks are added to the system color electric fields are created due to eq. H13a|l and one 
can distinguish two different situations: those with non-vanishing and those with vanishing total color charge. The 
prototypes of these two configurations are an isolated quark and a qq configuration, respectively. In both systems a 
bag with K « 1 in its interior develops which is stabilized by the vacuum pressure B. In the former case there exists 
only a monopole term and Gauss's law can be solved for the electric displacement Z?" exactly. Due to the radial 
symmetry the fields are perpendicular to the bag surface and is a smooth function of the radial distance r. Using 
radial symmetry the electric energy of the quark can be calculated 



^ - 47rr2 ^-1-2 7'^'^ «(r) " 8^ 7 r^^^ir) ' ^^^^ 

In our numerical realization the diverging self energy is regulated due to the finite quark width w{x). But the energy 
diverges also in the long range limit as soon as K,{r) vanishes more rapidly than pOj . In this case the vacuum 
pressure cannot balance the energy and the bag radius diverges as well. In the qq case the field lines start and end 
at the quark and the anti-quark. Thus they can arrange to be completely parallel to the bag surface and the electric 
field is a smooth function in space. The electric energy can be expressed as 

El!=\ I d\K{f)E^-E^ . (22) 

As K(r) vanishes rapidly outside the bag the energy density is localized within the bag and the electric energy stays 
finite. This geometry already can only be solved numerically. Solutions for the transverse profile of qq fields for large 
quark separations R with axial symmetry are given in [4^ . To illustrate the second scenario we will first present a 
qualitative picture based on a simple bag-like model, where the bag is a cylindrical tube with axial radius p and with 
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a sharp boundary, i.e. k = 1 inside and n = K^ac = outside of the bag, respectively. Afterwards the quantitative 
analysis will be based on the CDM where k becomes a smooth function of r. In the bag model the electric displacement 
vanishes exactly outside of the tube. Further we will assume that the electric field is homogeneous and constant 
inside the bag. The strength of the electric field inside is given by eq. H13a() as Z)" = E'^ = gsq"" /{irp^) and the 
total energy in a central slice of the tube with thickness A€ is i?tot = ^'ci + E^oh with E^i = ^olC^ / {tt p^) and 
i?voi = Birp^Ai. Minimizing the energy with respect to the radius p yields 

Pt = 9lCf/{2n^B) , (23a) 



Eo - El, + E% = j2glCfB M tM , (23b) 



El = K,y = ^glCfB/2M (23c) 

for the minimizing radius po and the corresponding energy Ef). Here we have introduced the concept of the string 
tension r. The identification of the string tension with the integrated energy density in the central slice of the string 
is only valid for strings with constant width, i.e. for infinite large quark separations R. As we will see, the radius 
increases for finite values of R. Therefore in the central slice the volume energy increases while the electric energy 
decreases. A more suited definition for the string tension in this case is r = , where i?tot is the total energy of the 
string including the end caps. At the stable point po the electric and the volume energy balance each other exactly 
(cf. eq. (|23c(l . For increasing values of B the dielectric vacuum compresses the electric flux into thinner flux tubes 
while for increasing electric flux, i.e. for increasing gs, the string radius grows (cf. eq^ 23a|) '). Note that the string 
tension scales linearly with the coupling gs which is typical for all bag models [43, 5C | . 

With a typical string tension r = 980MeV/fm taken from meson spectroscopy 22 > 113 and a tranverse radius 
of Po = 0.35 fm taken from lattice calculations we get a bag constant B^/'^ w 315 MeV and a strong coupling 
constant C^Us = Cf^ g1/{AiT) « 0.15. Though the value of the bag constant varies over a wide range in the literature 
from = 145 MeV in the MIT bag model Q to B^/^ = 241 MeV taken from QCD sum rule analysis [HEl, 

the value found here is rather high. On the other hand, the value of the coupling constant is rather small as 
compared to that obtained on the lattice Cpcts = 0.3 |53l | or in meson spectroscopy ranging from Cpcts = 0.3 — 0.5 
[32, 54]. Note however, that we have presented a qualitative discussion here based on a simple bag model with fixed 
and sharp boundary. We will study in section IIII Bl how the string tension r and the string radius p depend on the 
model parameters when the confinement field a and the electric fields are calculated according to the equations of 
motions (|16|) . For a long flux tube, when the quarks are located at the x-axis at x = ±a;o and separated by a large 
distance, one may assume axial symmetry for the geometry of the string. In the central plane between the quarks at 
a; = the fields can be described by uij) = a{p) and 0° = CaX. The constant Ca is determined by Gauss's law. In 
this case the dielectric displacement has the simple form 

D'^ip) = -Ca^i[p)e^- (24) 

The electric field points along the flux tube axis and its profile is proportional to the profile of k{p) = K,{a{p)). We 
shall see in sec. llVl if and for which quark separations R this asymptotic behavior is reached. The simple picture of 
sharp boundaries discussed above is reproduced if the dielectric function k{p) is a step function. In this limit the 
penetration depth of the electric field into the non-perturbative vacuum is zero. For any smooth profile of the 
dielectric function the penetration depth stays finite and non-zero but the electric fields are still screened. 



A. Generic features of CDM 



We proceed in studying in detail both the geometrical structure and the energetic content of a qq flux tube of flnite 
length. Our major interest is concentrated on the transverse profile of the energy density of this object as well as on 
the scaling of the total energy with the qq separation. These quantities will be compared later on to lattice results 
and experimental data. In this section we use a parameter set which we will call later PS-I and which is given in 
tab. IIVI below. The corresponding potential U{a) is shown in fig. (solid line). 

It should be noted that the numerical realization is not restricted to the qq geometry. The qq configuration possesses 
an axial symmetry, and the equations of motion Ijltill can be reduced to two dimensions. In the limit of an infinite 
long qq string, the problem even can be reduced to one radial dimension. But already the three-quark system does 
not have this symmetry and therefore we need a three-dimensional algorithm, as presented in appendix lA II 

To start with, we show the electric displacement of such a fiux tube in fig.^ The qq pair is located on the x-axis 
at X = ±0.6 fm. It is nicely seen that the field lines emerge from the quark as a source for the 3-component of the 
field and end on the anti-quark as a sink. Between the particles the field is nearly homogeneous and drops off at an 
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axial distance of about 0.3 fm. At the positions of the particles one sees the expected Coulomb-like field for point-like 
particle. In the same figure we show also the contour lines of the dielectric function K{a) for k — {0.2, . . . , 0.8} starting 
from the outside of the bag. Outside the bag where k < 0.2 the electric displacement nearly vanishes. 

It should be noted that the 8-component of this very special configuration (namely a rf pair) has exactly the same 
geometrical shape but is reduced by qf/Qr = ^/VS- Of course the electric fields for a different color/anti-color (e.g. 
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FIG. 4: The arrows represent the electric displacement & of a rr flux tube. The dots at a; = ±0.6 fm represent the (anti-) 
quark position. The solid lines are the equidistant contour lines of the dielectric function at «: = (0.2, . . . , 0.8) starting from 
outside. The string has a radius of p « 0.3 fm. The electric 8-field of this configuration is & = {q^ /q^)& . 

bb) pair are different as the blue quark has no 3-component of the charge (see fig. ^ . The two color configurations 
are connected by a color rotation V defined in eqs. and© which leaves the energy (or action) density invariant. 
The energy density is shown in fig. |5l At the quark positions the Coulomb peaks develop, but we do not show them 
here to emphasize the structure of the flux tube between the particles. The energy is distributed over a well defined 
region outside of which the energy vanishes. 

The profiles of the underlying fields D^, the confinement field a and the dielectric function k at x = of a 1 fm long 
qq string are shown in fig. El On the string axis k and a reach their maximal and minimal value, respectively, and both 
functions approach to their vacuum values Kvac and cTvac, respectively, for p — \/y'^ + > 0.5 fm. The confinement 
field does not drop down to cr = and therefore the fields are not completely in the perturbative phase. The region 
where k deviates substantially from its vacuum expectation value (p < 0.3 fm) defines the bag and consequently the 
electric flux falls off very fast at the bag boundary. 

In the same figure on the right we show the corresponding profile of the energy density. Here we have decomposed 
the total energy H14I) into an electric, a volume and a surface term according to eq. (|15|l . The surface energy (dash- 
dotted line) is maximal on the bag boundary at p « 0.3 fm, justifying the expression surface energy. The volume 
energy (dotted line) is distributed more homogeneously within the bag whereas the electric energy (dashed line) falls 
off quite rapidly. 

The energy density in the central plane between the charges at a; = is of special interest as it has been analyzed 
on the lattice in SU(2) '2^ and also in the framework of the dual color superconductor and in the Gaussian 

Stochastic Model 10, 56, _57j. Because of the symmetry all quantities depend only on the distance p from the string 
axis and we may reduce our analysis of the geometry of the flux tube to the shape of this profile function. We follow 
the reasoning of p9l | and compare the energy profile to both a dipole and Gaussian-like parameterization 

/d(p) = Nd{p^+pl)-' (25a) 
fg{p) = A^gexph In 2 (p/pg)"] , (25b) 

where n in Ij25b|l is a parameter giving the steepness of the profile. For n — 2 fg{p) is a Gaussian. For small quark 
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FIG. 5: Energy density (integrand of eq. HSU ) of the same qq string as in fig.|l] A well defined flux tube is stretched between 
the particles. The scale is chosen to show the flux tube and therefore the strong Coulomb peaks at the quark positions are not 
seen. 
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FIG. 6: The profile of the fields (left) and the corresponding energies (right) of a 1 fm string. The electric field is scaled to 
its central value on the axis and the scalar field a to its vacuum value. 



separations R perturbative QCD predicts Coulomb-like fields with a charaeteristic dipole behavior (|25a|l whereas at 
large quark separations the field should fall off much faster and the profile should be described better by a generalized 
Gaussian Ij25b|) with a width at half maximum pg. We show the profile of the energy density with the parameter 
set PS-I in fig. |3 together with the best fit of both parameterizations. One can see that the CDM results can be 
nicely described with a dipole shape for small separations {R ~ 0.4 fm) and with a Gaussian shape {n — 2) for large 
separations (i? ~ 1.0 fm). It should be noted that the Gaussian shape is only a qualitative guess for the profile. The 
profile might fall off even faster than described by a Gaussian. To quantify this we will extract below the parameter 
n from the fit of (j25b|) to the profile. A value of n = 2 indicates a pure Gaussian shape and a value n > 2 is connected 
to a sharper bounded bag. 

In addition to the shape of the flux tube, we will study below the dependence of the total energy (see eq. (|15|l ) as 
a function of the quark separation R, i.e. the qq potential Vqq{R). We compare the calculated CDM qq potential to 
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FIG. 7: The profile of tfie energy density togetlier witli tlie best fits of tiie parameterizations l|25|l fitted over tiie wiiole range 
sliown. Tfie profile is better described by a dipole form for small qq separations R (left panel) and by a Gaussian form for large 
R (right panel). In the left fig. we have omitted the dipole peak at the axis to better show the characteristic dipole tail. 



the Cornell potential 

Vc{R) = 2CfEo - ^ + tR. (26) 

This potential has been observed on the lattice in SU(2) and SU(3) 58] and has been used successfully in meson 
spectroscopy for heavy quarkonia [siL l32L Is^ . It contains the three parameters Eq, a and t which will be determined 
by a fit to the CDM results. For short distances R one might expect that perturbative one-gluon exchange results 
of QCD are dominant and thus Vc{R) shows the characteristic Coulomb-like y/R potential. To take into account 
non-perturbative effects a is an effective coupling constant. It has been shown |4J,E3 ^^^^ there is a 1/R correction 
in the gg-potential due to quantum-mechanical vibrations of strings around the classical solution. It is beyond the 
scope of the present paper to predict the corrections to the i?-dependence within our model. A possible correction 
of the Coulomb-type is included in the effective coupling constant a. For large R the linear term in eq. (|26|l takes 
over. It is a consequence of the formation of long linear flux tubes. The string tension t expresses the strength of 
confinement. The constant term is due to the self energy of the quarks. For point-like particles this should diverge 
but in our numerical realization we have regulated it by the non-zero quark width tq. The potential will vanish as 
i? — > 0, because the equal and opposite charge distributions of the quark and the anti-quark start to overlap and to 
cancel each other. Therefore a Cornell fit is only meaningful for quark separations R substantially larger than vq. 
For convenience we have separated from the Coulomb and the constant term the factor = 1/3. The agreement 
between our calculations and the Cornell parameterization in eq. ()26|l will be demonstrated in the next section. 



B. Variation of the CDM parameters 



In the remainder of this section we will discuss the influence of the parameters B, nig, cr vac and on both the 
energy profile and the total energy. For this purpose we vary all parameters one by one in a certain range while 
keeping the others at a fixed value given in tab. |nj We calculate the profile for a i? = 1.0 fm long string as well as 
the potential Vqq[R) for varying R. From the generalized Gaussian fit to the energy profile we extract the width pg 
and the steepness parameter n. From the Cornell fit to the CDM qq potential we extract the string tension r and the 
effective coupling constant a. The results are collected in tab. ITll 

We also study, if the equality between the electric part and the volume part of the string tension is still valid in 
our model. Note that this equality holds in the simple flux tube model expressed in eq. ()23c|) . For that purpose we 
also perform a Cornell fit to the different energy fractions Ec\, E^oi and i?sur separately. From these fits we extract 
the different parts of the string tension Td, TvoI and Tgur- We present the ratios /3 — Toi/tvo1 and 7 = Tgur/Ttot also in 
tab.m We find for all parameter combinations a remarkable equality between Tei and TvoI and the ratio /3 is 1 within 
a few percent. The surface part of the string tension amounts to roughly 7 = 0.2 except in the case when B ^ and 
all other model parameters kept fixed (see first panel in tab.^. In all simulations shown in this paper we have used 
Kvac = 10""*. All string quantities approach its limiting values very rapidly and do not change any more for value 
Kvac < 10~^ (see also appendix I A 2|l . 
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B^/^ [MeV] 




/3 


7 


Pa M 


n 





632 


1.07 


0.28 


0.44 


2.8 


60 


633 


1.07 


0.28 


0.44 


2.8 


120 


642 


1.07 


0.27 


0.43 


2.7 


180 


675 


1.06 


0.24 


0.40 


2.5 


240 


755 


1.06 


0.19 


0.35 


2.3 


nig [MeV] 




/3 


7 


pg [fm] 


n 


1000 


755 


1.06 


0.19 


0.35 


2.3 


1200 


831 


1.03 


0.19 


0.36 


2.6 


1400 


905 


1.00 


0.20 


0.36 


3.0 


1600 


982 


1.00 


0.21 


0.36 


3.4 


1800 


1053 


1.00 


0.22 


0.36 


4.0 


Cvac [fm"^] 




/3 


7 


Pg [fm] 


n 


1.01 


755 


1.06 


0.19 


0.35 


2.3 


1.26 


925 


1.06 


0.21 


0.35 


2.4 


1.51 


1110 


1.06 


0.23 


0.34 


2.4 


1.75 


1306 


1.06 


0.23 


0.33 


2.3 


2.00 


1511 


1.06 


0.24 


0.32 


2.3 


gs 




/3 


7 


Pg [fm] 


n 


0.2 


100 


0.94 


0.20 


0.21 


1.9 


0.5 


223 


0.97 


0.20 


0.24 


2.0 


1.0 


418 


1.06 


0.19 


0.29 


2.2 


1.5 


591 


1.06 


0.19 


0.32 


2.3 


2.0 


755 


1.06 


0.19 


0.35 


2.3 


2.5 


911 


1.05 


0.18 


0.37 


2.4 


3.0 


1062 


1.05 


0.18 


0.39 


2.4 


3.5 


1212 


1.04 


0.18 


0.41 


2.3 


4.0 


1356 


1.04 


0.18 


0.42 


2.3 


5.0 


1637 


1.04 


0.17 


0.43 


2.2 



TABLE II; The string tension r, the string width pg and the steepness parameter n for varying model parameters. The ratios /3 
and 7 are defined in the text. We vary the parameters one by one and keep the other fixed at B^^^ = 240 MeV, rUg = 1000 MeV, 
o"vac = 1.01 fm~^ and Qs = 2, respectively. Throughout this work we use Kvac = 10"*. 



a. Bag constant B: We start with a variation of the scalar potential U (see eq. {TJ). We vary B from B = 
to the maximal value allowed by eq. © B^^^ = 240 MeV (first panel in tab At the lower value U{a) has two 
degenerate vacua and at the higher value it has only an inflection point at cr = (see dashed lines in fig. From 
the simple bag model, eq. (|23|) . we would expect the string tension to increase quadratically with B^/^ and the string 
width to decrease with B^^^^ for increasing values B. Instead the string constant t grows only weakly over the whole 
range and is non-zero a,t B — 0. Note that for B = the profile of the volume energy U{p) is non-zero everywhere 
and a string of finite width is formed. The string radius shows the expected tendency though the decrease starts only 
for B^^^ > 120 MeV. The profile is steeper than a Gaussian for small B {n = 2.8) and approaches a Gaussian for 
larger values. 

b. Glueball mass nig: When we fix the bag constant and vary only the glueball mass nig (second panel in tab. lll|l 
the string tension increases roughly linearly with rUg. The string width stays constant over the whole tested range 
nig — (1000 . . . 1800) MeV in tab. m Simultanously the steepness n of the profile increases strongly up to n = 4. 
In the limit rUg — > oo the local maximum U{a) rizes infinetely and the confinement field a is restricted to be either 
(T = or (T = (Tvac • In this case we expect to reproduce the bag with sharp boundaries discussed in eq. . Of course 
this is not reached in the CDM but we observe the tendency in the rise of n. We see that the string tension is more 
affected by the mass mg than by the bag constant B and that vice versa the string radius is nearly independent on 
nig but decreases with B. Note that the glueball mass has no analog in the bag model, where the shape of the flux 
tube is assumed to be rectangular. 
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c. Vacuum value a vac- Qualitatively the string quantities should depend in the same way on the vacuum value 
CTvac as on nig. If we increase CTvac from its minimal allowed value given by eq. ^ to higher values a local maximum in 
the scalar potential U{a) develops (see dash-dotted curves in fig. El- In this case the stiffness rUg, i.e. the curvature 
of the potential at Uvac is unchanged. The increasing surface energy related to the larger gradient in the confinement 
field cr reduces the string radius pg. In the tested range the radius decreases from pg — 0.35 fm to pg = 0.32 fm (see 
third panel in tab.^l- Though the decrease in the string radius is not very large, the string tension doubles its value 
from 755 MeV/fm to 1511 MeV/fm. The steepness of the energy profile is rather independent from (Tvac- 

d. Coupling constant Qs'- The couping gs determines directly the strength of the electric flux in the string. From 
eq. H23|l we expect the string tension to rize linearly with gs and the string radius to increase proportional to ^/gl■ 
Indeed we find a nearly perfect linear behavior of t (4th panel in tab. ^ and also the string width has the expected 
qualitative behavior. The linear dependence between the string tension and the coupling gs^fCp is valid in general 
for all bag-like models [H^l- For the coupling g^ we also look at the qq potential at small quark separations i?, 
i.e. we extract the parameter a in the Coulomb term of the Cornell potential which is dominant for i? < 0.2 fm. For 
small R the bag has almost spherical shape and the qq interaction can be treated perturbatively. In fig. |51we show 
the Coulomb parameter for varying gs- The electric part ad (circles) is perfectly described by the perturbative result 

2 

Q!oi = (dashed line). The surface contribution asur (triangles) does not follow the quadratic behavior and seems 
to grow more linearly and rather slowly. It thus acts as a non-perturbative correction to the perturbatively expected 

2 

Coulomb interaction a — |^ + asur- The volume energy increases linearly with R and thus the corresponding Coulomb 
parameter is compatible with zero and not included in the figure. 




FIG. 8: The Coulomb coupling parameter cx extracted from the Cornell fit for different gs- The electric part (dots) is in 
perfect agreement with the one-gluon exchange result (line). The surface part (triangles) contributes only a small amount to 
the Coulomb potential and therefore the effective a (squares) rises quadratically with gs. 
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TABLE III: The reaction of string tension r, the string width pg and the steepness n on increasing model parameters. 



We end this section in summarizing the reaction of the string tension r, the string width pg and the steepness of 
the profile on the parameter variation in tab. IIIII 
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IV. qq STRINGS 
A. CDM and lattice results 

After having analyzed the influence of the model parameters on the shape and the energy of qq strings, we compare 
the CDM results to lattice calculations. We first describe qq strings to fix the parameters and later on turn to qqq 
systems in section Ivl 

In detailed studies of long flux tubes were made in SU(2) lattice gauge theory. The authors compared the 
profile of the action and the energy density of strings of varying length. We will tune the parameters of our model 
to reproduce the phenomenological value of the string tension r and to desribe the energy profile of a 1 fm long qq 
string as good as possible within our model. For the profile a measure for the agreement of our calculations and 
lattice results is the quantity (Ae)^ = X]i[^(Pi)~£iat(/Oj)]^; where e\at{pi) is the energy density given at discrete radial 
distances p,; (see fig. 24 in To find the optimal set of parameters we scan a wide range of our parameter space 

and minimize Ae with the constraint to reproduce the string tension r — 980MeV/fm. 

Doing this we follow three different strategies. In the first we treat all model parameters from tab. IIIII as free 
parameters, determined from the minimization of Ae only. In this case we do not lean on the heuristic interpretation 
of B as the bag constant and nig as the glueball mass. This results in parameter set PS-I in tab. II VI In the second we 
take B and rUg as the bag constant and the glueball mass, respectively. The former is chosen to be B^/^ — 240 MeV 
as found in QCD sum rules |^ [s^l and also in another bag- like analysis js^ ■ The latter is restricted to be at 
rUg = (1500— 1700) MeV as given by SU(3)-lattice results for the scalar glueball mass 0,|S2j. This leaves only gs and 
Cvac as free parameters and rUg in the small range mentioned to find the minimum of Ae and we get parameter set 
PS-II in tab. lIVI In the last variation we restrict additionally the coupling constant to gs = 3.3 in order to reproduce 
not only the string tension but also the Coulomb parameter a — 0.29 of the Cornell potential as calculated on the 
lattice |53| . The last free parameter cr vac is varied in this case to get the right string tension and we obtain parameter 
set PS-III. Of course the deviation of the CDM energy profile to the lattice result increases from parameter set PS-I 
to PS-III. The three different parameter sets therefore express a quantitative variation of the model predictions. With 
the given parameters the scalar potential U{a) changes as shown in fig. |^ Within PS-I (solid curve) there is only 
a negligible relative maximum in the potential, and it vanishes exactly for PS-III (dash-dotted line). In the second 
parameter set a pronounced relative maximum develops (dashed line). We note that the values of the bag constant 
B are much larger than the value B^^'^ = 145 MeV chosen in the original MIT bag model T| and also in a previous 
CDM analysis (2^]. Recall that in the string profile is much broader and inconsistent with lattice results. The 
vacuum value for the dielectric constant is Kvac = 10"'' for all parameter sets. We show in appendix lA 21 that the 
physical quantities of the string do not change anymore for smaller values of Kvac- 



No. 


B 


nig 


^ vac 


9s 




I 


(260 MeV)* 


1000 MeV 


1.29 fm"^ 


2.0 


10-* 


II 


(240 MeV)* 


1500 MeV 


1.13 fm"^ 


1.8 


10-* 


III 


(240 MeV)* 


1700 MeV 


0.59 fm-^ 


3.3 


10-* 



TABLE IV: CDM parameter sets used in the description of qq strings and qqq baryons. 

In fig. 1101 we show the result of the three different fitting procedures for the 1 fm long qq string. The profile for 
PS-I runs smoothly through the lattice points showing the high quality of the fit. The profile of PS-II has roughly 
the same half maximum width = 0.32 fm but a steeper profile due to the higher glueball mass. In the last case 
for PS-III the flux tube is much broader with a half maximum width pi/2 = 0.45 fm. Due to the large coupling gs for 
PS-III we cannot push the string radius to smaller sizes while keeping the string tension at the prescribed value. 

To analyze further the different parameter sets, we decompose the total energy according to eq. (|15|l into the 
different energy parts. The result is shown in the profiles in fig. 1111 In all cases the electric energy builds up most 
part of the total energy. Note that this is not in contradiction to the previous result that the electric and the volume 
part of the string tension are of the same magnitude. The string tension is by definition r = and not the energy 
of the central slice of the string. This would be the same only for strings with constant width, which is not the case 
for finite quark separations. 

For PS-I (left panel) the volume energy (dash-dotted lines) never exceeds the bag constant, whereas for PS-III 
(right panel) it has a wide range inside the bag, where it is constant and equal to B. For PS-II (central panel) we see 
that on the string axis, the volume energy comes close to the local minimum at u = 0, which is seen in the central 
dip. The surface energy is strongly pronounced for PS-II and the two peaks are clearly separated for PS-III. In the 
last case the interior of the bag is therefore much more pronounced than in the other two parameter sets. 
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FIG. 9: The scalar potential U{a) for the three parameter sets given in tab. IIVI 
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FIG. 10: The energy density profile of a 1 fm string. The profiles for PS-I and PS-II have a half maximum width of p ~ 0.35 fm 
although the slope is steeper for PS-II. PS-III results in a broader string. Lattice results taken from j29l| . 



The underlying fields D°-{p), a{p) and the dielectric constant K,{p) for this 1 fm string are shown in fig. El For all 
parameter sets the scalar field never reaches the perturbative situation cr = within the string. According to that n 
is never exactly equal to 1. For PS-I and PS-II it has only a value of k = 0.75 and k = 0.85 respectively. Only for 
PS-III the perturbative value k < 1 is reached. Of course also the fields are broader as compared to the other two 
parameter sets. 

The dielectric constant has roughly the same shape as the electric field D°'{p). This indicates already that 
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FIG. 11: The decomposition of the energy profile into the different energy components for K=\ fm. 
according to eq. (|24|l the string picture should be vahd for i? = 1 fm. One can expect that for increasing qq distances R 
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FIG. 12: The scaled profiles of D", a and of a 1 fm string. The profiles of k and Z?" have roughly the same shape as 
expected from eq. I|24|l for long flux tubes, k reaches 1 only for PS-III and is smaller than 1 for PS-I and PS-II. 

the profiles evolve to some stable shape and that the asymptotic relation eq. (|24(l between electric field and dielectric 
constant should become increasingly accurate. To study this issue we show the profile of the total energy density for 
different quark separations R in fig.E) For small R we see the strong Coulomb peak. But for separations i? > 1.2 fm 
the string profile does not change strongly anymore and the asymptotic profile is nearly reached. 

This can be seen as well in fig. El where we have plotted the ratio of the normalized profiles of and k in the 
range \y\ < 0.7 fm, where the string is located. According to eq. H24I) this should be constant equal to unity for large 
R. Indeed this ratio becomes increasingly flat and equal to one for increasing R. 

To complete the discussion of the energy density profile we compare the full half maximum width A — 2pg extracted 
from the Gaussian fit to the results obtained in lattice SU(2) calculations 29] for various string lengths R. Here we 
have fixed the steepness parameter to n = 2 to be consistent with the analysis in |22j. In fig. one sees that the 
calculated string widths are compatible with the lattice data for parameter sets PS-I and PS-II. Again the string 
width is overestimated for PS-III due to the large coupling g^. First the width increases rapidly for small separations 
R < 0.8 fm. For the largest quark separation i? = 2fm the width has a value A = 0.78 fm and A = 0.68 fm for 
parameter sets PS-I and PS-II, respectively, and the width seems to saturate. The width for PS-III is A = 0.97 fm 
and still increases slightly. We note that qualitatively this behavior is consistent with the lattice string picture [g^Ig^ . 
where a logarithmic increase of the string width with the quark separation R is predicted. 
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FIG. 13: The profile of the energy density for different qq separations 7? given in the figure. The asymptotic shape is nearly 
reached for R > 1.2 fm. For small R the Coulomb peak is seen. 




FIG. 14: The ratio of the electric field to the dielectric constant k, each normalized to its central value. The asymptotic 
relation in eq. 12411 is nearly reached for qq separations R > 1.2 fm. 



B. CDM and the Dual Color Superconductor 



In this part of the discussion of qq strings we make a connection to the Dual Color Superconductor model (DCS) 
known also as the dual Ginzburg-Landau model. In this model color flux tubes between a quark and an anti-quark 
are formed by the interaction to a scalar field carrying a magnetic charge. Those Abrikosov-Nielsen-Olesen vortices 

..65. ,66(1 are stabilized by a circular magnetic current 

Jmag V X ^DCS (27) 

flowing around the string axis. It is speculated that such a magnetic current is produced by the condensation of 
magnetic monopoles which can be constructed in non-Abelian gauge theories |3 13 • The DCS flux tube has a 
characteristic profile of the electric field -Edcs — -£'DCs(p)ex- To compare our results with the DCS model we first 
note that by construction of the electromagnetic field tensor F'^'^'" in eq. I)ld|) the curl of the electric field E vanishes 
identically for static configurations. However, by a simple redefinition of the field tensor F^'^'"' — > QP-f,a _ ^(fj^^/^i'.a 
one may define a quantity 



(28) 
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FIG. 15: The Gaussian width of the flux tube as a function of quark separation. For parameter sets PS-II the width saturates 
at A « 0.7 fm. For PS-1 and PS-III the width of the profile still increases slightly for separations up to J? = 2fm. Lattice data 
are taken from |29|. 



where Q^^'^''^ — L^oif3tiiyQtj.iy,a ^^iq dual field tensor to G^'^'"^ . This quantity exactly behaves as a magnetic current and 
the spatial part is in our standard notation J^^^g = V x D". In the absence of magnetic fields there is no magnetic 
charge density, i.e. p^j^g = ^ ' = 0- o^'' model the magnetic current is connected to the electric displacement 
D"" whereas it is connected to the electric field E-qcs- Thus l5° and i?DSC are confined in the two models and we 
compare i5° with S^cg . In (GTj, J3£] a detailed study of the electric profile of a qq configuration was made within 
lattice gauge theory and it was found that the profile can be described very well by the DCS model. In addition the 
authors found a magnetic current Jmag fulfilling eq. H27|) . In fig. ^| we show the lattice results of the electric profile 
(triangles) compared to those obtained within our calculations. We find only a qualitative agreement, where again 
PS-I and PS-II reproduces the electric profile better than PS-III. 

The spatial part of the magnetic current constructed in eq. H28(l in the central plane between the two particles at 
a; = is displayed in fig. El (left). In a cylindrical basis it has only an azimuthal component. One sees the circulating 
structure of the current which gives this fiux tube the name vortex in the DCS model. We show our results for the 
profile of the magnetic current with the different parameter sets and compare it to lattice data [g^ in the right panel 
of figure El The maximal values are shifted to larger values of p compared to lattice data. Parameter set PS-III 
develops a pronounced current only on the surface of the string. In identifying the magnetic current in our model we 
can make a link between our model and the model of the dual color superconductor. It should be noticed that the 
CDM model is formulated on the basis of the gauge potentials A° which can be directly related to the gluon fields of 
QCD, whereas the DCS model is formulated in the dual gauge potentials. In this sense the model gauge fields in the 
CDM can be more easily interpreted as the QCD gluon fields than in the dual superconductor model. 



C. The qq potential 



We can quantify the analysis of a string further by showing the total energy E'tot of the string as a function of 
the qq distance R in fig. 1181 together with the electric, the volume, and the surface part of the energy. All parts of 
the energy show a linear rise with the qq separation for R > 0.2 fm. Due to the self energy of the nearly point-like 
particles the electric energy is larger than the other two energy contributions. This is most obvious for PS-III where 
the coupling constant gs is largest. We have performed a Cornell fit according to eq. (|26|l to the total energy i?tot 
and to the single energy parts separately. The extracted parameters Eq , a and t and the corresponding values for the 
electric, volume and surface energy contributions are listed in the first three columns in tab. ]V\ It should be noted. 
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FIG. 16; Comparison of the profile of the electric displacement D within CDM and the profile of the electric field E found on 
the lattice. 
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FIG. 17: The circular magnetic current >nag = V x Z?" (leftl and the scaled profile of it (right). For PS-III the current is locate 
at the surface of the bag. Lattice results taken from Wt\ . 



that the self energy i?o depends on the width of the particles. The given values correspond to a Gaussian width of 
ro = 0.02 fm which can still be resolved on the used computational grid. 

The string tension is r « 980MeV/fm as the model parameters were fitted to this value. The effective strong 
coupling ap = C^a isolated from the Cornell fit to the qq potential ranges from 0.12 with PS-II to 0.30 with PS-III. 
A Cornell fit obtained from meson spectroscopy gives an effective Coulomb coupling ranging from ap = 0.25 |54| to 
Op = 0.5 (32]. The latter estimate included both the charm and the bottom quark mesons. In the MIT bag model 
fw^ and in the CDM hadronization study an effective coupling ap = 3.0, and ap = 2.6 were used, respectively. 
Together with the small value of the bag constant used in these works the string width amounts to p = 1.9 fm and 
p — 1.5 fm, respectively, which is large as compared to the value of p = 0.35 fm obtained on the lattice .29i] . In a 
further bag model analysis the authors extracted a value ap — 0.38 for the coupling constant. The coupling 
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isolated from a Cornell fit to lattice data prefers a value of ap = 0.3 [ssL l61*|. The CDM values for the coupling in 
tab. IVl underestimates this value for ap for PS-I and PS-II. Parameter set PS-III was designed in order to reproduce 
the Cornell potential, so the agreement is optimal here. 
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FIG. 18: The qq potential E{R) for the different parameter sets together with the different fractions of the energy and the 
Cornell fit to the total energy. 
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PST 


PS-II 
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PS-II 
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PS-III 


CD [MoV] 


388 


294 


757 


961 


763 


2202 


2.5 


2.6 


2.9 


ao 


0.18 


0.12 


0.30 


0.41 


0.31 


0.87 


2.3 


2.7 


2.9 


T [MeV/fm] 


979 


982 


980 


1548 


1502 


1586 


1.6 


1.5 


1.6 


eg [MeV] 


297 


235 


726 


833 


679 


2163 


2.8 


2.9 


3.0 




0.12 


0.09 


0.29 


0.34 


0.28 


0.86 


2.7 


3.1 


3.0 


Tel [MeV/fm] 


404 


390 


424 


640 


585 


695 


1.6 


1.5 


1.6 


Tvoi [MeV/fm] 


381 


380 


435 


610 


588 


717 


1.6 


1.5 


1.7 


e'S' [MeV] 


108 


78 


47 


151 


108 


63 


1.4 


1.4 


1.4 


sur 


0.07 


0.05 


0.03 


0.10 


0.07 


0.04 


1.4 


1.4 


1.4 


T,ur [MeV/fm] 


195 


212 


122 


299 


328 


174 


1.5 


1.5 


1.4 



TABLE V: The parameters e_D — CdEq, an = Cna and r extracted from the Cornell fit to the total energy as well as to the 
different energy contributions separately. The first block of three columns belongs to quarks in the fundamental representation 
{D = F), the second to quarks in the adjoint representation {D = A) and the last shows the ratio of both. 



After having shown the CDM results for the qq potential, we would like to comment on the so called Casimir 
scaling hypothesis 0, [7^ [7^ 0, [t^- According to this hypothesis the qq potential should scale with the 
quadratic Casimir operator Cu of the static sources q and q in the representation D . For small quark separations 
R the potential is dominated by the perturbative Coulomb term which already scales with Cp). For larger distances 
this cannot be deduced from perturbation theory and must be shown numerically on the lattice. This has been done 
for various representations D in ,.72,. jd. ItgI Ft^ . For the adjoint representation for example one should expect the 
string tension to scale with {Ga/Cf)tf = 2.25tf. This can hold if at all only for intermediate quark distances, 
as the adjoint potential should saturate for separations larger than some critical distance Rc (Ta. l77l |. Numerical 
values for the string tension tjj have been given in 72] and in [tsI ]. An unambiguous confirmation for the scaling 
hypothesis has not been seen. The deviations from scaling where found to be (2-5)% and (10-15)% [l^ for various 
representations. Interestingly, all results lie consistently under the value predicted by Casimir scaling. For the ratio 
of the adjoint string tension to the fundamental one a value ta/tf = 1.97 ± 0.01 ± 0.12 was calculated in |7i|, with 
the respective statistical and systematic errors. This corresponds to a deviation from Casimir scaling by (12 ± 6)%. 

In the CDM we can simulate adjoint quarks. We simply assign the adjoint charge to the three quarks, i.e. the 
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sum of the charges of a quark and an anti-quark. The numerical values are given in tab. IVII Note that one can 
construct only 6 charged adjoint quarks. The other two members of the octet are uncharged. We have already seen 
that the Coulomb term in our potential scales with the square of the coupling constant gs or, rephrased in terms of 
the Casimir operator, scales with the Casimir Operator. The string tension, however, was found to scale with gs and 

y^C^, respectively. Therefore the adjoint string tension should be ta — \J C^'/C'^ Tp fa l-lrp, where we have used 

the Abehan Casimir values = 1/3 and = 1- 

We have performed the simulation for the adjoint string and isolated the Cornell parameters ca — CaEq and 
aA = Cao: and ta- The results of these values for the different parameter sets are given in tab. [V] (2nd block of 
three columns) together with the corresponding ratios (3rd block). The short range Cornell parameters of the electric 
energy have approximately the expected ratio of 3. Also the adjoint string tension as well as the different parts of 

the energy scale with the expected value of ^ C^/C^ = 1.7 which has to be compared to the SU(3) scaling value 

Ca/Cf = 2.25 and the empirical value ta/tf = 1.97 W^. In SU(2) lattice calculations a clear deviation from 
Casimir scaling was detected. We conclude that Casimir scaling is not seen in our model which is in accordance with 
all bag models 0j E0| ■ Casimir scaling is also not seen within the Dual Color Superconductor model, when the 
superconductor is on the border between first and second type [SIZ^, i.e. where lattice results agree with this model 
the best. We note that it is also qualitatively in line with the tendency of lattice calculations to underestimate the 
adjoint string tension systematically. 

We note in this context that the deviation is an intrinsic feature of the CDM, which is due to the larger fluxtube 
cross section in the present of charges in higher representation. If in the future lattice data will show that Casimir 
scaling would hold for the gg-potential, this would be a strong argument against all bag-like models in favor of models 
with explicit scaling 10] . In the meantime we regard those type of models as competitive. 



color q q 

rg 1 

rb 1/2 ^3/2 

gb -1/2 73/2 

gf -1 

bf -1/2 -V3/2 

bg 1/2 -73/2 



TABLE VI: Color charge of color sources in the adjoint reprontation. 



V. qqq BARYONS 

Now that we have fixed the model parameter to reproduce lattice results fo r qq string s, we go on to describe 
baryon-like qqq configurations. It has been discussed for a long time |3a. 1371 133. l39l l57l l79 | . whether the flux tubes 
stretching between the quarks will connect the particles pairwise with a qq string, indicating 2-particle interactions, 
or whether they are connected via a central point giving rise to a real 3-particle force. In the former case the geometry 
of the system will show a triangular or A-like shape, whereas in the latter it will have a Y-like shape. The situation 
is depicted in figure [T^ Although our model is formulated in the Abelian approximation we must not naively assume 
that the flux tubes are a simple superposition of three qq flux tubes. The non- linear interactions with the dielectric 
medium might deform the flux tube to show a Y-like geometry. This question has been studied on the level of the 
potential within lattice SU(3) 39, 5^ [s^l and on the level of the fields [s^l. Depending on the shape the 

potential will scale characteristically with the 2-particle distance R. To parameterize it in the spirit of the Cornell 
potential eq. we decompose the potential into a constant term due to the quark self energy, a Coulomb-like short 
distance term and a confining linear term. The constant term scales with the number of particles in the system and 
the short range term scales with the sum over the two particle Coulomb interactions. The Coulomb interaction is 
accompanied by the same color factor ^eaPj^aTyt^a^pr — as in eq. (j20bp . 

The confining term scales with the total length of the flux tubes spanned between the quarks, which is different 
in the two geometries (see fig. I19II . In the case of the Y-geometry a qq-like string is connected to each of the three 
quarks, meeting at a central point. In the case of three quarks sitting on the corners of an equilateral triangle with 
length R, the length of each of the flux tubes is equal to R/VS. This yields eq. H29a|l below. 
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FIG. 19: A-shape (left) and Y-shape (right) of the qqq configuration. In the Y-geometry each of the three quarks is connected 
to the central point with the same string as in the qq case. In the A-geometry the quarks are connected pairwise by a modified 
flux tube. 



In the A-geometry two flux tubes are connected to each quark, so there is a reduced electric flux in each qq string 
compared to that of the qq string. Due to the symmetry of color exchange the total energy per unit length stored in 
the two strings must be of equal size. To deduce this relative strength of the modified flux tube we apply a limiting 
procedure. Think of two of the quarks approaching each other, e.g. the r-quark and the g-quark. As they come 
together we end up with a quark-diquark system but with two qq strings lying on top of each other. As the diquark 
behaves exactly like an anti-quark the linear part scales with the known string tension t and we conclude that the 
qq flux tube has a string tension reduced by a factor of two. This factor is also obtained in [s^. Note that we have 
assumed here, that the flux in each qq string is independent of the third quark position. The qqq potential may now 
be parameterized as: 

W = iCfE, ^^Y. + - -'^l (29a) 



3CfSo-3^f- + V3Ti? 
z K 



Here are the positions of the particles and fc is the position of the flux tube junction in the Y-case. The parameters 
Eq, a and r are the ones extracted from the fit to the qq string in the previous section and listed in tab. The 
respective last equality in eq. (|29|l holds for the quarks sitting on the corners of an equilateral triangle. The effective 
string tension ta = (3/2) r is smaller than ry = VSt making it the preferable configuration although it is only an 
effect of 14%. 

This picture is surely oversimplified, since the flux tubes are not strings with zero transverse extent. The single 
flux tubes will overlap due to their finite width p and the true fields will lead to a smearing of the two extreme 
configurations at least for quark distances only slightly greater than the flux tube width. Both the finite size of the 
flux tubes and the tininess of the overall effect makes it hard to decide on the lattice, which might be the more 
appropriate description. In references [sl, |5^ a A-like scaling of the potential was found suggesting an effective 
2^article interaction, whereas in [3^l83 | the baryonic string tension was better described by Y. In another calculation 
135] neither ansatz give a proper description of the potential for all quark separation. Instead the authors stated that 
the potential is of the A-type for small separations R and of F-type for large R. 

In the following we first show the fields for the three-quark system with the symmetry shown in figure 1191 as 
calculated in CDM. Then we will compare the qqq potential obtained within CDM to those constructed in eq. 129f) 
using the Cornell parameters Eq , a and r obtained in section IIVI Finally we will perform a new fit of our Cornell 
parameters to the qqq potential. We would like to stress here, that in a renormalization group derivation of the 
dielectric model [T^lldis^ an additional colorless vector field appears. This vector field is relevant in the calculations 
of systems with non- vanishing baryon-density. However in the spirit of the present work, we are interested in the pure 
glue qqq-type fiux tubes with fixed external charges, where the neglect of this additional term is justified. 

In figure QUI we show the color electric field of a configuration, where the rightmost, uppermost and lowermost 
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quarks have color blue (&), red (r) and green {g), respectively. The pairwise quark distance is equal to i? = 1.7 fm. In 
the left (right) part of the figure the electric 3-field (8-field) D^^^ is shown. As the b quark has no 3-charge component 
the 3-field connects only the r and the g quark. But as can be seen clearly, due to the existence of the b quark, the 
confinement field a is reduced in the whole region between the three quarks and the electric 3-flux is bent towards 
the b quark. The electric flux of the 8-field connects all three quarks. The r and the g quark are equal sources of 
this field and the flux ends on the b quark. Again the flux is deformed compared to the qq flux tube and is pushed 
towards the center of the qqq configuration. We have included in the same picture the contour lines for the dielectric 
function K{a). One can see that the maximal value of k (neglecting the peak positions at the quarks) shows a Y-type 
shape, and falls off to zero towards the sides of the triangle within a range of 0.5 fm. 




FIG. 20: The color fields of a qqq baryon. The red/green/blue quark is the upper/lower/rightmost one. The flux (left) 
connects only the red and the green quark, whereas the flux (right) connects all three quarks. The flux is bent into the 
center of the conflguration. Solid lines are the equipotential values of k. 



As we have argued in section ^ the electric field shown in figure 1201 is not an invariant quantity under the global 
color symmetry. The strength and the direction of the fields depend on the specific arrangement of the colors in the 
qqq state. We have included the picture of the electric field only to give an insight to the underlying mechanism 
of flux tube formation. The relevant quantity again is the energy density of the system. As in the qq case we 
split the total energy into electric, volume and surface terms H15II . In figure we show the energy density of the 
different types together with the total energy density obtained with parameter set PS-I for a qqq state with quark 
separation R = f .7fm. A clear difference in the geometries is seen for the different energy fraction. The electric energy 
distribution (upper left panel) shows pairwise electric flux tubes between the three quarks, bent into the center of the 
baryon. The geometry is nearly A-shaped due to the dip at the center. Going along the x-axis from the center into 
the negative direction, there is a well in the energy density of about 95 MeV/fm"^, corresponding to 20% of the central 
value. This is qualitatively the same for all parameter sets, but the energy barrier is larger in PS-II (f 20 MeV/fm'^) 
and smaller in PS-III (OOMeV/fm^). 

In contrast to the electric part of the energy the volume part has a clear Y-shaped structure (upper right panel) 
which is the same as for the dielectric function k((t) shown in figure EOl The surface part (lower left panel) is only 
relevant where a varies spatially. This is true on the edges of the triangle and one sees therefore a pure A-like 
distribution. The absolute value of the surface energy density in the flux tubes is small compared to the values in the 
electric and the volume fluxes. 

The sum of all three energies results in the total energy density shown in the lower right panel. The picture shows 
qualitatively the same structure as the one for the electric energy. But the barrier of the total energy is only 10% of 
the central value and the competing structures of the electric energy and the volume energy are smeared out. 

We will study now whether this picture of superimposed structures can be found in the baryonic potential when the 
size of the qqq state is varied. The result is shown in figure together with the Y and the A-parameterization of the 
Cornell potential (|29|l . The absolute values of the total energy agree very nicely with the A-parameterization which is 
in line with the results in |38l l53| . However, it is the string tension Tqqq that distinguishes the two parameterizations 
and not the absolute value -Btot- If one takes a closer look at the potential, one sees a slightly larger slope than 
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FIG. 21: Energy density of the qqq configuration. Black dots show particle positions. Wfiite voids at particles are spikes in the 
energy density due to the particles Coulomb peaks which are cut out. Upper left: eei, upper right: Cvoi, lower left: esur, lower 
right: etot- Lines are equidistant in energy density with the values given in the figures. 
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FIG. 22: The total energy of the qqq state compared to the Y and the A-forms of the Cornell parameterization. 
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expected from the A-picture. To analyze the baryon potential with respect to Tqqq we make a fit of the form 

3 a 

Vqqq = iCp Eq - -Cp + Tg^^i? (30) 

to the CDM results and extract the string tension Tqqq from the fit. Here Eq and a are free fit parameters as well. As 
in the qq case we make this fit to the total energy (I15a|l as well as to the electric, the volume and the surface parts 
of the energy (|15bp - (|15d|l to see if the different contributions show the same behavior as the energy distributions in 
figure EH 
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PS-II 


PS-III 
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1544 
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0.72 
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278 
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St 
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-0.016 


0.01 



TABLE VII: The qqq string tension of the baryon Cornell potential and the deviation St from the expected Y {St = 1) and 
A-picture {St = 0). 



The result of the fits is summarized in table IVIll Here we show only the values for the qqq string tension as both 
the constant term and the Coulomb term scale in the same way in the A and the Y-geometry. Together with Tqqq 
we show the deviation St ~ {Tqqq — ta)/(ty — ta) from the expected values in both pictures. A value St = (1) 
indicates perfect agreement with the A (Y) picture whereas St « 0.5 means that neither the Y nor the A is realized 
but a transition between both geometries. 

First we confirm the visual impression that the electric and the surface part of the string tension t°^^ and r^^^ are 
well described by the A-picture. All deviations from the A parameterization St are close to or even undershoot zero. 
In the surface part of PS-I this can be explained by the respective energy density (lower left panel in fig.[2J- There 
most part of the energy is located outside the region defined by the triangle of the quarks. The opposite is true for 
the volume part r™^. Here the deviations are close to one, indicating a better description by the Y-geometry. In the 
qqq case the equality between the electric and the surface part of the string tension is not true anymore. For PS-III 
Tqqq IS larger than Tqqq by 20% with the same tendency in the other two parameter sets. 

The total string tension Tqqq as a superposition of the three parts takes on a value in between the A and the 
Y-value. All parameter sets have a slight preference for the A-picture with parameters St between 0.21 (PS-III) and 
0.41 (PS-II). 

Another useful check for the scaling is to look at the qqq-potential for arbitrary geometry of the ggg-triangle. With 
the parameterizations for the Y and the A-picture one would expect a universal behavior of the long range part of 
the potential with the total string length L. For the Y-geometry we have Ly = J^i l^i — '^cl and for the A-geometry 
La = J2i<j ~ For triangles with no angle exceeding 120° the minimal string length L\ can be calculated |8J| 
with 

1/2 

^Y= U V(r;-f,)2 + 2x/35A I , (31) 



where 5'a is the area of the triangle spanned by the three quarks. We show in fig. 1231 the qqq potential for arbitrary 
quark positions. We have selected only those triangles, where all angles are smaller than 120° and where the 2-particle 
distance \ri — fjl > 0.4 fm. The visual impression agrees with the previous result, that the ggg-potential is better 
described by the A-geometry. The potential plotted as a function of Ly scatters more strongly in the energy than in 
the A-case where it follows a rather straight line. Note that this result is different to that obtained in [s^l . 

In summary the qqq configuration shows two competing geometries visible in the energy distribution and in the 
potential. The scalar field has a Y-like structure while the electric energy follows a A-like distribution. In the former 
case the bag energy is minimized by minimizing the volume of the three-quark bag. However, the flux tubes have a 
finite transverse extent with no sharp boundaries. Therefore the electric flux is not restricted to the Y-shaped flux 
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tube but may evolve in a region outside the scalar Y-volume. The superposition of both the electric and the scalar 
energy distributions leads to a potential right in-between the two ansatze, which is qualitatively the same result as 
obtained in a lattice calculation ,8^ . 



VI. CONCLUSIONS AND DISCUSSION 



In this work we have shown that the Chromo Dielectric Model is able to describe confinement of overall color 
neutral systems by the formation of color flux tubes in a perfect dielectric vacuum. We studied the dependence of the 
string tension and the string profile on the model parameters. The parameters of the model were chosen carefully to 
reproduce the profile of qq strings as well as the qq potential obtained in lattice calculations and from heavy-meson 
spectroscopy. Three different parameter sets were found that describe the profile of the qq flux tube obtained on 
the lattice optimally under given constraints. All parameter sets reproduce the phenomenological value of the string 
tension r — 980MeV/fm. Parameter set PS-I additionally describes both the width as well as the shape of the flux 
tube rather well. In PS-II we fixed the bag constant and the glueball mass parameters to values from the literature 
and found that only the width of the profile can be reproduced. The shape develops a steeper profile than the 
lattice profile. In a final set we have fixed as well the coupling constant to reproduce also the Coulomb parameter 
of the Cornell potential. With this parameter set PS-III the resulting string width is larger by 20% than the value 
obtained on the lattice. The profiles of the strings reach their asymptotic shapes for quark separations larger than 
R ~ 1.2 fm. However, the width does not saturate at a constant value but increases slowly, which is in line with the 
lattice expectation, where the string width increases logarithmically. It would be desirable if there were SU(3) lattice 
data available of the same accuracy as obtained within SU(2) lattice calculations. Then one could compare the CDM 
results to a somewhat more realistic theory. 

We have also made a comparison of the qq flux tubes obtained in CDM to the calculations within the Dual Color 
Superconductor. We were able to extract a magnetic current in the CDM showing the same vortex-like behavior as 
in the DCS. Both the proflle of the electric field and the magnetic current were in qualitative agreement with the 
profiles calculated within the Dual Color Superconductor model and also calculated on the lattice. 

Finally we have shown that the qq string tension does not obey the Casimir scaling hypothesis. The scaling of the 
adjoint string tension with ^/CaJCf and not with Ca/Cf is a general feature of all bag-like models. Casimir scaling 
has not been proven unambiguously on the lattice and therefore does not exclude our model. 

For qqq configurations we have calculated the color electric fields, the color invariant energy density and the qqq 
potential to discuss the geometric form of those baryonic states. We have found two competing pictures in the electric 
and the scalar part of the system. The scalar energy distribution is confined into a clear Y-like form whereas the 
electric energy distribution is of the A-type form. In the total energy the overall values are in good agreement with 
the A-type of the qqq Cornell parameterization but the string tension Tqqq shows the same competing behavior as 
the energy distributions: A-like in the electric sector and Y-like in the volume sector of the energy. The total string 
tension thus has neither the Y nor the A-like value but lies rather in-between the two pictures. 
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Out of the three parameter sets PS-I to PS-III the first one describes best all lattice data of qq flux tubes shown 
in the present work, although the values of the model parameters differ from that found in the literature. It should 
be noted, that the identification of nig with the glueball mass and of B with the gluon condensate obtained in QCD 
sum rules is motivated from heuristic arguments only. Thus the results obtained with PS-I give a good agreement 
to lattice data and the other two serve for a comparison if the space of model parameters is restricted to the values 
given above. The Chromo Dielectric Model gives therefore an adequate mechanism of confinement, which can be 
understood dynamically. 

It would be interesting to see the influence of the color electric fields within a self-consistent treatment including also 
the quarks dynamically. Those calculations were done in |l3j| including a direct interaction between the confinement 
field and the quark field only. In our model the interaction with the quark field is only indirect via the quark-gluon 
interaction — g g'ip^ ^t"" A^^ij} and the interaction of the gluons with the confinement field. Here is the quark field 
obeying the Dirac equation. 

A natural extension of our calculations is to include also time-dependent fields according to eq. ((T^ . Those 
studies were done in where only the confinement field was treated dynamically. Within those time-dependent 
calculations one could study the hadronization out of a gas of colored quarks and gluon fields which might be produced 
in a relativistic heavy-ion collision. Hadronization is up to now treated as an instantaneous process at the freeze-out 
temperature. Within our model one could test this assumption as a time-resolved process. Of course this is numerical 
challenging due to the complexity of the equations but might be feasible within nowadays computing power. 

APPENDIX A: NUMERICS 
1. The FAS- Algorithm 

To solve the set of equations we define a rectangular box with fixed volume V — y. Ly x = and 
discretize the equations on a Cartesian grid Vth with N — it' nodes and n — 2™ — 128. The grid spacing is therefore 
h = L/n. The particle positions are not restricted to this grid, but can be placed arbitrary. To avoid spurious 
oscillations in observable quantities when changing the quark position, we assign a spatial width to the quarks given 
by tq. For the grids used in this work with i = 3fm we have a grid spacing h w 0.02 fm and this determines the 
width ro = 0.02 fm. 

We use the Full Approximation Storage (FAS) multigrid algorithm described in [s^, IH, 1131 . The FAS algorithm 
is especially suited for our purposes because both the needed memory resources and the computing time scale only 
with the number of grid points N which is mandatory for those large systems. It is like any other multigrid algorithm 
an iterative solution technique that improves on an initial guess for the solution step by step. Conventional relaxation 
methods cease to converge for grids with increasing n (decreasing h). This means that the convergence get worse if 
one wants to decrease the discretization error. To be more precise it is the low frequency part of the system which 
does not converge anymore. The grid introduces an infrared cutoff due to the finite dimensions L of the box and 
an ultraviolet cutoff due to the finite number of grid points n (known as Brillouin zone in solid state physics) . The 
Fourier spectrum of the solution is built up of modes with wave vector ki — (jijL . . . ii: / L . . . fcmax — wk jV). Modes 
of the error made in approximating the true solution with wave vector ki and i> njl (oscillating modes) are damped 
out fast, but modes with i < n/2 (smooth modes) do not die out due to the locality of the dilference operator. 

To cure this problem the algorithm introduces temporary coarser grids fie with ne — {n/2, n/4, . . . , 2~^ri, . . . , 2) 
nodes per dimension. On those grids the spectrum is reduced to modes with maximal wave vector kf^^^^ = 2~^/cmax 
and smooth modes on the fine grid become oscillating modes on the coarse grid. Further relaxation of the solution 
on the coarse grid allows to solve for those modes and in turn to get an estimate for the error made on the fine grid. 

Thus the FAS algorithm consist of three building blocks: (i) the relaxation method on each grid, (ii) the transport 
of the approximation from a fine grid to the next coarser grid and (iii) the transport back from the coarse to the fine 
grid. The recursion starting from the finest grid flo = flh down to the coarsest flm-i and back is called a V-cycle. In 
contrast to linear multigrid algorithms the FAS algorithm allows also for non-linear differential equations. 

Both the electric potentials 0° and the confinement field a are solved within the same algorithm at once. As a 
smoothing method we use /i = 4 Gauss-Seidel relaxations with red-black ordering for the smoothing update on each 
grid both on the downward and the upward stroke of the V-cycle. For the discretized version of eq. Ijl6b|) we have 
to include a Newton-Raphson approximation in order to cope with the non-linearity. The transfer of all discretized 
field quantities from coarse to fine grids and back is performed with linear interpolation and full weighting reduction 
operators respectively. To improve the convergence of the algorithm we have to find a good initial guess to start the 
V-cycle on the finest grid JIq = ^h- To this end we first find a solution on the next coarser grid fii which is much 
less expensive in the computational costs. In general we find a solution on every grid Qi with an initial guess found 
on the grid fle+i- 
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The computational costs of the multigrid V-cycle in d space dimensions can be estimated as foUows |8g • The 7 = 3 
independent fields (u, 0") and the corresponding residuals defined below in eqs. I|A1|I are stored on a grid with n'^ 
nodes. Thus the minimal amount of memory needed is given by Mi — mijn'^ fioat value units with nii — 2. On 
each of the n'^ nodes the discretized equations for the 7 fields are solved locally 2/i times and the 7 residuals are 
calculated once. If we estimate the amount of work for the transfer between the fine and the coarse grid with another 
2n'^ operations, the minimal number of computational operations on each grid is given by M2 = {2fi + 3)'yn'^ = m2^n'^ 
with 7712 = 11. We may call the minimal amount of memory Mi and the minimal number of numerical operations 
M2 one memory unit and one working unit respectively. The equivalent numbers on a coarse grid fi^ are M^^^ — 

TOi/27(?i/2^)''. An upper limit for the numerical costs summed over all coarse grids can be estimated by a geometrical 
series X^fco -^^1/2 ~ ™i/27"-''(l ~ 2*^)^^ = |toi/27?^'^ for d = 3. The algorithm is surprisingly more efficient in higher 
dimensions d. In d = 3 space dimensions the amount of memory needed is therefore only times larger than that 
for traditional relaxation methods and the computational work for a whole V-cycle sums up to |77i2 times that for 
one relaxation sweep. One sees that the total computational costs scale only linearly with the number ^n'^ unknowns. 

We impose Dirichlet boundary conditions on the confinement field a — dvac and von-Neumann boundary conditions 
on the electric potentials, i.e. = where n is a normal vector to the boundary. This forbids any electric flux out 
of the computational domain and is in line with the observed fact, that the electric flux is compelled into the flux 
tube. 

We compute after each V-cycle the residuals 

r% = V-Ka)V0'^)-Kg.p" (Ala) 
r„ = - [/'(a) + -E . (Alb) 

with the approximations 0° and a and the corresponding energies given in eq. (|15|l . We observe a rapid decrease of 
\r^ \ and \r„ \ and simultaneously a fast convergence to the asymptotic values of the energy. We show as an example r„ 
and i?tot as a function of the computational work in figure !^ for the multigrid algorithm (solid symbols) and for pure 
Gauss-Seidel relaxation (open symbols). For both methods we have started with the same initial guess obtained on the 
next coarser grid. On the left panel one sees, that the rate of convergence is slowing down for Gauss-Seidel relaxation 




FIG. 24: The residual (left) and the total energy (right) as a function of the number of V-cycles in the FAS algorithm. 
Gauss-Seidel relaxation are the open symbols, the FAS results are given with solid symbols. Note the logarithmic scale for the 
computational work. 



when going from small grids (n = 32) to larger grids (n = 128) whereas it is independent of the discretization within 
the FAS algorithm. On the right panel we see that the total energy converges very rapidly to the asymptotic one in 
the FAS algorithm. In practice we stop the algorithm if the total energy stays constant within a given error of 10~^. 
This is reached usually after three to five V-cycles on the largest grid. The time needed for the solution is reduced 
by nearly 2 orders of magnitude as compared to ordinary relaxation. On a 2.2 GHz Pentium 4 machine a typical 
configuration is solved on a — 128"^) grid within 5-6 minutes. 
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2. The vacuum value of the dielectric constant 

In the spirit of the model one should set Kvac — 0. However the Poisson eq. H16a() is ill-defined in the perturbative 
phase in this limit. In 47] it was shown that in order to get a consistent solution for the flux tube, one also must 
have K'(crvac) = 0- Therefore we have parameterized the dielectric constant according to eqs. (|10|l and with a 
finite but small value for Kvac- With this form also K'(crvac) for Kvac — > 0- As we are dealing with a boundary 
value problem we could not start right away with K'(crvac) = as in this limit a = cr vac everywhere is a solution of 
(jl6b|l which will be found by the algorithm. Therefore we have chosen a rather large value Kvac = 10~^ when finding 
an initial guess on coarser grids with n < 32. We then decrease the value of Kvac step by step on the finer grids until 
a prescribed minimal value is reached. We have analyzed the dependence of the string quantities on /^vac- 

In fig. EH 

we show that the string tension as well as the profile parameters do not depend on Kvac, once it is smaller than say 
Kvac — 10 . The numerics were still stable and fast for Kvac — lO"'^ which we have used throughout this work. Of 
course it is not possible to calculate the energy of a single-quark configuration, as the radius of the bag diverges and 
most of the electric energy is stored outside the bag. But we can estimate this energy in the same spirit as for the 
cylindrical flux tube with sharp boundaries in eq. and obtain Eg — ^[{2g'^CFB) / {tt"^ k^^^)]^^'^ . With the smallest 
values used in our model ^ 240 MeV and g,, = 2) we obtain Eq ~ 66 GeV which is large on a hadronic scale. 




FIG. 25: All parts of the string tension (left) and parameters defining the shape of the profile (right) saturate very fast for 
decreasing values of Kvac- In all results shown in this work, we have used a value Kvac = 10~*. 
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